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THE  NEGATIVE  BINOMIAL  DISTRIBUTION: 

COMPUTATION  OF  THE  MEDIAN  AND  THE  MEAN  ABSOLUTE  DEVIATION 


Definitions 


For  n>  0,  0  <  p  <  1,  and  q  =  1  -  p,  the  distribution  of  the  discrete 
variable  x,  having  frequency  function 


f(x;n,p)  - 


x-tn-  lx 

\  n  x 
x  p  q 


-n 


N  n  .  .  x  ^  _ 

x  )  p  (-q)  ,  x  -  0,1,2 . 


i 8  called  the  negative  blnomia 1  distribution .  It  also  is  called  Pascal 1 s  distri¬ 
bution  when  n  is  a  positive  integer,  and  is  called  the  geometric  distribution 
when  n  ■  1 . 

If  n  is  a  positive  integer  there  are  two  well-known  instances  of  this 

distribution.  In  a  sequence  of  Bernoulli  trials  with  probability  p  of  success, 

t  h 

f(x;n,p)  1  s  the  probability  that  the  n  success  will  occur  on  the  trial 

numbered  (nbx) ;  that  is,  it  is  the  probability  that  exactly  x  failures  will 

precede  the  n^  success.  Also,  f(x;n,p)  is  the  frequency  function  of  the 

distribution  of  the  sum  of  n  random,  independent  variables,  each  of  which  has 

the  geometric  distribution  with  frequency  function  f(x;l,p);  that  is,  f(x;n,p) 

t  h 

is  the  frequency  function  of  the  n  convolution  of  the  geometric  distribution 
with  Itself,  which  will  be  evident  from  the  generating  function.  DOC 
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Mean,  Variance,  and  Hlgner  Moments 
Mean ,  p  =  n  q/  p 
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Variance  from  the  mean,  a  “  n  q/ p 


"\ 
J  j 


uu 


SEP  2  3  1966 


C 

These  values  can  be  obtained  by  direct  summation  or  from  the  generating  function 
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Fv*>  "  )  f(x;nkp)  s*  ■  *  pri(l-qt)  ° 

x-0  ' 

Higher  moments  can  be  computed  from  the  formula 

i 

f  ("Tl)(f)k  . 

x-k 

obtained  by  taking  the  kth  derivative  of  F^s)  and  putting  s  -  1.  Such  moments 
are  needed  in  fitting  polynomials  in  exponential  smoothing.  Thus,  if 


pm  ’),  \  <T)  tk  • 

00 

the  n^  exponentially-weighted  average  of  the  values  of  the  polynomial  for 
t  <c  0  is 


00 
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x  f(x;n,a) 


However,  if  the  prlynomial  is  written  in  the  form 
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It  it  for  this  reason  that  the  second  font  of  the  polynomial  is  the  preferred 
form. 

Median 

In  general,  the  equation 

m 

f<x;n,p)  -  1/2  , 

x-0 

has  no  Integral  solution  m.  However,  just  as  in  the  positive  binomial  distri¬ 
bution,  the  partial  sum  can  be  vritten  in  terms  of  the  Incomplete  Beta- function, 
which  then  can  be  used  as  the  definition  of  the  partial  sum  for  non-integral 
values  of  m.  In  this  way  we  can  find  a  value  of  m,  not  necessarily  an  Integer, 
that  satisfies  the  equation.  This  value  of  m  will  be  called  the  median. 

The  formula  for  the  positive  binomial  distribution  is 

m 

n  ^  pX  qn  X  *  I  (n-m,mH)  ■  1  -  I  (ntH,n-m),  (1) 

x  )  y  H  H  P 

x-0 


where 


is  the  Incomplete  Beta-function.  The  formula  appears  in  many  books  and  can  be 

proved  easily  by  integration  by  parts. 

The  corresponding  formula  for  the  negative  binomial  distribution  is 


m 
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This  formula  la  not  readily  available.  It  is  stated,  but  not  prominently,  in 
the  Introduction  to  Pearson's  Tables  of  the  Incomplete  Beta- function  and 
Pearson  gives  a  proof  (provided  by  Fieller)  in  Biometrika,  Vol.  XXV,  pp.  160-161. 
A  simpler  proof  is  the  following:  Integrating  by  parts, 
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Hence 


m 

t  /  B  n  V-1  /  r&m\  m-k  k 

Ip  (n.mfl)  “  p  ^  V  k  J  P  q  <3) 

k“0 

By  induction  on  m  it  is  easy  to  show  that 


k-0  k-0 


Formula  (2)  is  obtained  from  (3)  and  (4). 

Although  formula  (2)  has  a  meaning  only  when  m  is  a  non-negative  Integer, 
the  integrals  in  the  Incomplete  Beta- function  exist  for  non-integral  values  of 
m.  We  define  the  median  of  the  negative  binomial  distribution  to  be  the  solution 
ra  of  the  equation 


I  (n,mfl)  -  1/2  (5) 
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A  unique  solution  m  0  exists,  provided  pn<C  1/2. 
Mean  Absolute  Deviation 

The  mean  absolute  deviation  from  the  median  is 


Let 
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x-m  f(x;n,p) 


[m]  ■  integral  part  of  m 

■  largest  integer  that  does  not  exceed  m. 
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Other  forms  for  A  can  be  obtained  from 

[m]  W 
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Two  of  these  are 


A 


(|i-ra) 


1-2  Ip(n,  [m]  +1) 


n 
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and' 


(m  -®) 


l"2lp(n,  [m]+l)J 


xp(n»  [m]  +2)  -  Ip(n,  [m]  +  1) 


(6) 


The  latter  form  is  easy  to  use,  since 


m 


+  1  m  +  1  <  ra  +2; 


that  is,  the  two  arguments  involved  are  the  two  Integers  between  which  we  inter 
polate  in  finding  the  solution  of  (5).  Thus,  to  find  A  t  enter  the  tables 
of  the  Incomplete  Beta-function  and  record  the  values 


I  (n,  I  ml  +1)  and  I  (n, 
p  l  J  p 


for  which 


yn, 


m  -HI)  <  1/2  and  I  (n, 

P 


+  2)  >  1/2 


when  the  second  argument  has  integal  values.  Interpolate  to  find  m  for  which 

I  (n,m+-l)  °  1/2  and  then  substitute  in  (6). 

P 

If  m  is  an  integer, 

A  (m+l)[\,/  ,  ,  1  _  o  ( n-1  mHl 
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If  m  is  not  an  integer  and  we  use  linear  interpolation  between  integral  values 
to  find  it,  then 


A  -  (  p  -  m 


+  it-K 


p(m-  [mi ) 


1  -  2  I  (n,  [m]  +1) 


A  quantity  of  interest  in  inventory  problems  is  the  expected  amount  bv 
which  x  exceeds  a  given  value  k,  that  is,  the  expected  back-order 


00 


B  -  )  (x-k)  f(x;n,p)  . 

L _ i 

x“k 


By  the  same  arguments  used  to  find  A  we  find 


B  -  (p  -k)  [  1-1  (n,kfl) 
L  P 


+  0±k) 

P 


1  (n.kf 2)-I  (n,W-l) 
P  P 


for  the  negative  binomial  distribution. 

Examples 

The  values  of  A  and  o  listed  in  the  table  below  were  computed  primarily 
to  test  the  hypothesis  that 


A  ■  k  o  , 


where  k  is  0.75  approximately.  For  p  °  0.1  it  is  necessary  to  use  the  formula 


I  (a ,b)  -  1  -  I  (b,a) 
P  <1 
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For  example 

•  xo. 

9 ( 6 , 1 )  -  0.5314  and 

I0.9^7, 

and  I0>1(1, 

7)  “ 

0.5217. 

2 

n 

2 

a 

0.1 

1 

9 

9.5 

2 

18 

13.4 

3 

27 

16.4 

4 

36 

19.0 

5 

45 

21.2 

0.5 

1 

1 

0.71 

2 

2 

2.00 

3 

3 

2.45 

4 

4 

2.83 

5 

5 

3.16 

0.9 

10 

1.11 

1.11 

20 

2.22 

1.57 

30 

3.33 

1.92 

40 

4.44 

2.22 

50 

5.56 

2.48 

“  0.4783;  from  which  Iq  ^(1,6)  ■  0.4686 


m 

H-m 

A/  a 

5.6 

3.4 

0.69 

14.5 

3.5 

0.75 

23.4 

3.6 

0.76 

32.4 

3.6 

0.77 

41.4 

3.6 

0.77 

0 

1 

0.71 

1 

1 

0.75 

2 

1 

0.76 

3 

1 

0.77 

4 

1 

0.78 

0.43 

0.68 

0.88 

1.51 

0.71 

0.85 

2.63 

0.70 

0.83 

3.74 

0,70 

0.82 

4.85 

0.71 

0.81 
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